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Abstract 



We present an observation and rate measurement of the decay D° —* K + it~it 
produced in 9 fb _1 of e + e~ collisions near the Y(4f>) resonance. The signal 
is inconsistent with an upward fluctuation of the background by 4.9 standard 
deviations. We measured the rate of D° — ► K + it~it normalized to the rate 
of ~W -> K + tt~tt° to be 0.0043±g;ggiJ (stat) ± 0.0007 (syst). This decay can 
be produced by doubly-Cabibbo-suppressed decays or by the D° evolving into 
a D° through mixing, followed by a Cabibbo-favored decay to K + ir~Tr°. We 
also found the CP asymmetry A = (9^22)% to be consistent with zero. 
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The transition of a D° into a D° through mixing provides a window through which we 
may observe the effects of non- Standard Model physics. Just as K° - K° and B° - B° 
mixing gave prescient information about the charm and top quarks before their discovery, 
observation of D° — D° mixing could imply evidence for new particles as massive as 100-1000 
TeV f . 

In this Letter we report the first observation of the "wrong sign" decay D° — > K + 7r tt° 
(consideration of charge-conjugate modes is implied throughout this Letter). The flavor of 
the initial D° was tagged by the sign of the slow pion, 7r s , from D* + — > -D°7r+. We measured 
the ratio R of wrong sign (WS) to right sign (RS) D° — > Kim decay rates, integrated 
over decay times. The WS decays can be produced by mixing of the initial D° to a D°, 
followed by Cabibbo-favored decay (CFD) to K + n^n , or by doubly- Cabibbo-suppressed 
decay (DCSD). 

The transition of a D° to a D° can proceed through real or virtual intermediate states, 
which we describe by the normalized amplitudes — iy and x, respectively 0. The Standard 
Model contribution to \x\ is suppressed by at least tan 2 9c ~ 0.05 due to weak couplings, 
however GIM cancellation || could further suppress \x\ by one to four orders of magnitude. 
While the Standard Model contributions are most likely below the present experimental 
sensitivities, many non-Standard Model processes could lead to \x\ > 0.01 ||. 

The interference between the mixing and DCSD amplitudes is influenced by a strong 
interaction phase difference between the CFD and DCSD amplitudes. We denote the ratio 
of DCSD and CFD amplitudes by — ^/R~£>e~ l5 . The leading minus sign is motivated by the 
relevant Cabibbo-Kobayashi-Maskawa matrix elements: Ti{V c dV* s /V cs V* d } < 0. An average 
over relevant three-body configurations of K + ts:~tx is implied in Rd and in the strong phase 
5[§. 

The effect of the strong phase is to yield measurable mixing amplitudes y' = ycosS — 
x sin 5 and x' = x cos 5 + y sin 5. Then, 



R = -±= = R D + \ R D y' + ~(x' 2 + y 

Data for this measurement were produced in e + e~ collisions within the CESR ring at 
center of mass energies near the T(4S) resonance. Data corresponding to 9 fb _1 of inte- 
grated luminosity were collected using the CLEO II. V upgrade of the CLEO II detector 
between February 1996 and February 1999. Reconstruction of displaced decay vertices from 
charmed particles was made possible by the improved resolution of the silicon vertex detector 
(SVX) 0, installed as a part of this upgrade. We utilized this improved resolution in previous 
searches for D° — D° mixing |8j and in measurements of charmed particle lifetimes PJlOU . 



Candidates for D° — > K + 7t~tt° were formed by combining good quality charged tracks 
detected in the drift chamber with 7r°'s formed from pairs of photons detected in the Csl 
crystal calorimeter. Photons from the central (end) region of the calorimeter with energies 
greater than 30 MeV (60 MeV) were considered. The tt° candidates were required to have 
momentum greater than 340 MeV/c, diphoton invariant mass within two standard deviations 
of the known ir° mass, good mass fit chi-squared, and at least one photon detected within 
the central region of the calorimeter. The tracks from the charged decay products of the D° 
candidate were required to form a vertex with confidence level of at least 0.0001. The track 
from the tt s candidate was refit with the constraint that it pass through the intersection of the 
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D° candidate direction and the CESR beam spot. This refit was required to have a confidence 
level of at least 0.0001. D* candidates with momentum less than 2.5 GeV/c were rejected. We 
separated signal from background using distributions of D° candidate mass, M = m(K7nr°), 
and energy released in the D* + decay, Q = m(K7CTT°TT s ) — m(A7T7r°) — m w . Charged kaon 
and pion daughters of the D° were required to have specific ionization consistent with their 
respective hypotheses. Combinations from RS decays in which both the charged kaon and 
pion were misidentified were removed by requiring the mass of the interchanged charged 
track hypothesis, m(irK7r ), to reconstruct at least four standard deviations away from the 
known D° mass. Similar kinematic vetoes were applied to m(KK7r°) and mlmrn ) in order 
to remove combinations with a single particle misidentification. 

In this analysis, systematic uncertainties were reduced by directly fitting for the scale 
factor S that relates the large number of RS events, Nr$, to the modest number of WS 
events, Nws- Nws = S ■ Nrs- Then, R = C ■ S, where the correction C can deviate from 
unity because the WS events can populate the Dalitz plot differently than the RS events do, 
and thus have a slightly different average efficiency. 

The scale factor S was determined using a binned maximum likelihood fit to the two- 
dimensional distribution of Q and M. The prominent and nearly background-free RS signal 
in the data was scaled by the factor S to provide the WS signal distribution for the fit. 
Background Kti7t 7t s combinations were broken down into three categories according to their 
expected distributions in Q and M: 1) RS D° — > K + tt~tc° daughters combined with an 
uncorrelated slow pion, 2) charmed particle decays other than correctly reconstructed RS 
D° — > K + 7t~7t°, and 3) products of e + e~ — > uu, dd, or ss events. Events from bb production 
were excluded by the D* candidate momentum cut. The contribution from RS D° — > 
K + 7t~7t° combined with an uncorrelated slow pion produces a peak in M, but is smooth in 
Q. While some backgrounds produce peaks in one variable, none of them produce peaks in 
both Q and M. Q-M distributions for the three backgrounds were taken from Monte Carlo 
simulations, however their normalizations were allowed to vary in the fit. We generated 
approximately 40 million e + e~ — > uu, dd ss and cc Monte Carlo events, corresponding to 
approximately eight times the integrated luminosity of the data, using the GEANT-based 



CLEO detector simulation 11 



The fit to the WS signal in the Q-M plane determined a scale factor S of 
0. 0043jlQ;ooio (stat) between the RS and WS signal peaks. This corresponds to a WS yield 
of 38 ± 9 events within the signal region of two standard deviations about the known Q and 
M values. The statistical significance of this signal was evaluated by fitting with the signal 
contribution constrained to zero and comparing the maximum log-likelihood value with that 
of the nominal fit. The signal was found to be inconsistent with an upward fluctuation of 
the background by 4.9 standard deviations. Projections of the WS signal and fit results in 
the two fit variables are shown in Fig. |l|. 

The correction factor C was estimated by fitting the WS Dalitz plot. This correction 
can differ from unity if the RS and WS Dalitz plots have different resonance structure, 
since the efficiency is not uniform across the Dalitz plot. Recently, CLEO performed a 
Dalitz analysis of the RS decay D° — > K + n~n° 0] and found a rich structure consisting of 
p(770)-, K*(892)+, A^(892)°, p(1700)", A (1430)°, A (1430)+, and A*(1680)+ resonances 
and non-resonant contributions. The dominant intermediate state in the RS decay is D° — > 
A+p(770)-, followed by W -> A*(892)+vr- and W -> A*(892)°7r°, which account for 
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roughly 79%, 16%, and 13% of the yield, respectively. The sum of these is greater than 
100% due to interference. 

We used an unbinned maximum likelihood fit to extract the relative contributions of the 
three major resonances from the WS data. Due to the limited statistics and relatively large 
background, only the amplitudes and phases of D° — > A"*(892) + 7r" and D° — > A*(892)°7r° 
relative to D° — > K + p(770) _ were varied in the fit. Relative amplitudes and phases of the 
minor contributions were fixed to the RS values found in [12| . 



The fit used the signal fraction from the Q-M fit and parameterizations of the efficiency 
and background distributions in the Dalitz plot variables as inputs. An analytic expression 
for the efficiency function was obtained by fitting a large sample of non-resonant signal 
D° — > K + 7T~7T° Monte Carlo events. The background function was estimated by fitting 
side-band regions in Q from the WS data. These regions contain contributions from RS 
D° — > K + 7t~7t° combined with an uncorrelated slow pion, e + e~ — > uu, dd, ss, and non- 
signal charm events. The measured RS squared amplitude JT2J, multiplied by the efficiency 
function, was used to describe the Dalitz plot of the RS D° — ► K + n~n° with uncorrelated 
tt s contribution. The combined e + e~ — ► uu, dd, or ss and non-signal charm fit function was 
taken to be a two-dimensional polynomial with coefficients determined from the Q side band 
fit. The relative contributions of these backgrounds were obtained from the Q-M fit result. 

The results of this fit were surprisingly consistent with those of the RS fit flI2| , leading to 
a correction of C — 1.00 ± 0.02 (stat). While the statistical uncertainties on the amplitudes 
and phases were large due to the low statistics and large backgrounds, the correction factor 
C was insensitive to these because of the relatively uniform efficiency. This resulted in a 
statistical uncertainty on C that was small compared with that for S. 

The total systematic uncertainty on R was estimated to be 17%. Contributions to this 
uncertainty are categorized in Tab. | as measurement errors on S and C from the Q-M and 
Dalitz plot fits, respectively. 

The important systematic uncertainties on S were due to uncertainties in the Monte Carlo 
simulation of the background Q-M distributions (14%), sensitivity to misidentification of 
charged D° daughter tracks (3%), and the statistics of the Monte Carlo sample (2.4%). These 
are summarized in Tab. |I[ The largest of these, due to the simulation of the background Q-M 
distributions used in the fit, was estimated by varying the Q-M side band regions used to 
constrain the background contributions. The sensitivity to misidentification of the daughter 
K + and 7r™ was studied by observing the variation of S with changes in the cuts applied to 
specific ionization and the kinematic cuts applied to m(irK7r°), m(KKir°), and m(7T7T7r ). 
The uncertainty due to the statistics of the Monte Carlo background Q-M distributions 
was estimated by performing a series of fits in which the contents of the bins were varied 
according to Poisson statistics. The uncertainty was taken to be the standard deviation of 
R based on 4000 variations. 

Uncertainties in the measurement of C come from the unknown minor resonance contri- 
butions to the Dalitz plot (8%), the Dalitz plot fit method (3.6%), the unknown background 
Dalitz plot shapes (3%), and the statistical uncertainty in the Dalitz plot fit (2%). These 
contributions are summarized in Tab. |I[ The amplitudes and phases of resonant and non- 
resonant components of the WS signal other than the dominant A + p(770)~, A*(892) + 7r~, 
and A*(892)°7T modes were fixed to the RS values in the fit. In order to explore the uncer- 
tainty of this assumption, we allowed these to vary in the fit. The systematic uncertainty in 
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the Dalitz plot fit method was determined by fitting the WS Dalitz plot under hypotheses 
that the signal was composed entirely of K + p(770)~, K*(892) + 7r~ , or K*(892)°7r° decays. 
C was found to differ from one by 2%, 1%, and 28%, respectively. The uncertainty was 
estimated by evaluating the consistency between the WS data and the pure K*(892)°tt° 
hypothesis, which produced the largest deviation from unity. This hypothesis was found to 
be inconsistent with the data by 7.7 standard deviations. The systematic uncertainty due 
to the Dalitz plot of the background was estimated using a series of fits with background 
shapes from side bands in Q and M, obtained from both the WS and RS data. When using 
side bands in M, the kinematics of the daughter tracks of D° candidates were scaled to force 
the allowed phase space to be similar to that of a true D° — > Ktttt decay. The statistical 
uncertainty on C from the Dalitz plot fit was included as a systematic uncertainty on R. 
The Dalitz plot fit results were checked by performing Q-M and Q fits to specific Dalitz plot 
subregions dominated by a single submode. The relative yields from these fits were com- 
pared and found to be in agreement with the efficiency-corrected squared amplitude from 
the Dalitz plot fit, integrated over the same subregion. The total systematic uncertainty on 
C was estimated to be 9.5%. 

By performing the analysis separately for D° and D° candidates, we measured the CP 
asymmetry of this decay, defined to be 
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R(D° -> K+ir~n°) - R(D U -> K 
R(D° -> K+tt-tt ) + R(D° -> K-tt+^ 



We observed an asymmetry consistent with zero: A = (9^.^)%. Due to cancellation of errors 
in this ratio, the systematic uncertainty in this measurement was negligible compared with 
its statistical uncertainty. 

In summary, we observed a signal for the decay D° — > K + tt~ , k q using 9 fb _1 data collected 
with the CLEO II. V detector. The signal is inconsistent with an upward fluctuation of the 
background by 4.9 standard deviations. This result is the first observation of the WS signal 
D° — ► K + 7t~7t°. Using fits to the Q-M and Dalitz plots, we measured the normalized WS 
rate and CP asymmetry to be R = 0.0043l{j;gftS (stat) ± 0.0007 (syst) and A = (9±1)%, 
respectively. 

To allow comparison with previous measurements of DCSD and mixing in the D° sys- 



tem |3,13 -15 1 we plot a band corresponding to this measurement of R in the Rd~v' plane, 
shown in Fig. The band depends on \x'\ and we show it for \x'\ = and \x'\ = 0.028, 
which correspond to the limits from our previous analysis of D° decay to K + n~ || if equal 
strong interaction phase differences are assumed for the decays ||. If we assume that there 
is no mixing, this measurement corresponds to Rp = (1.7 ± 0.4 (stat) ± 0.3 (syst)) • tan 4 Q C - 
We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Research Corporation, the Natural Sciences and 
Engineering Research Council of Canada and the Texas Advanced Research Program. 



7 



TABLE I. Systematic uncertainties in the R measurement. 



Measurement 


Source 


Uncertainty 


S 


MC background Q-M dist. 


14% 




K/ir separation 


3% 




MC statistics of Q-M dist. 


2.4% 


c 


Minor resonances 


8% 




Dalitz fit method 


3.6% 




Background Dalitz plot 


3% 




Dalitz fit stat. error 


2% 


R = C-S 




17% 
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FIG. 1. Results of the Q-M fit to the WS data, shown in projections onto a) Q and b) M 
within the signal region of two standard deviations about the known M and Q values, respectively. 
Contributions to the WS data sample come from the WS signal (WS), RS D° decays combined 
with an uncorrelated slow pion (RS), decay products of e + e~ — > uu, dd, or ss (uds), and decays 
from charmed particles, other than correctly reconstructed RS events (Charm). 
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interaction phase differences are assumed for the two modes. The bands indicate the region within 
one standard deviation of this measurement. 
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